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We study vacuum polarization of quantized massive scalar fields (j) in equi- 

' librium at black-hole temperature in Reissner-Nordstrom background. By 

o ■ 

' means of the Euclidean space Green's function we analytically derive the 

• renormalized expression < 0^ >h at the event horizon with the area 47rr^. 
^ ' It is confirmed that the polarization amplitude < cp'^ >h is free from any 

divergence due to the infinite red-shift effect. Our main purpose is to clarify 
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• the dependence of < cp"^ >h on field mass m in relation to the excitation 
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' mechanism. It is shown for small-mass fields with mr_|_ <C 1 how the ex- 
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citation of < 0^ >h caused by finite black-hole temperature is suppressed 
as m increases, and it is verified for very massive fields with mr-|- ^ 1 that 
< cf^ >H decreases in proportion to with the amplitude equal to the 
DeWitt-Schwinger approximation. In particular, we find a resonance behav- 



^ ' ior with a peak amplitude at mr_(_ ~ 0.38 in the field-mass dependence of 



vacuum polarization around nearly extreme (low-temperature) black holes. 
The difference between Scwarzschild and nearly extreme black holes is dis- 
cussed in terms of the mass spectrum of quantum fields dominant near the 
event horizon. 
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I. INTRODUCTION 



Quantum behaviors of matter fields in black hole spacetimes have been extensively stud- 
ied for understanding the various physical effects. In particular, the existence of a state of 
quantum fields in equilibrium at a finite temperature near the event horizon has attracted 
much attention, because it clearly represents the thermodynamic properties of stationary 
black holes. The problem of vacuum polarization for this Hartle-Hawking state may be 
described in terms of the Euclidean space Green's function Ge{x, x'), which is periodic with 
respect to the Euclidean time r = it. If one considers a quantized scalar field 0, the vacuum 
polarization < > can be calculated by using the equality 

< (j)\x) > = Re{lim Ge(x,x')} , (1) 

x' — >x 

in which the renormalised expression is derived through the method of point splitting. 

It is well-known that the black-hole temperature T defined as the inverse of the period 
of Ge( 

X J X ^ IS proportional to the surface gravity k on the event horizon as follows, 

T = k/271 . (2) 

(Throughout this paper we use units such that G = c = h = ks = !■) If the origin of 
the vacuum polarization < > is claimed to be purely induced by the finite black- 

hole temperature, the amplitude should decrease toward zero in the extreme black-hole 
limit K ^ 0. In fact, we can see this behavior of < 0^ > by applying the analytical 
approximation of the renormalized value obtained by Anderson, Hiscock and Samuel 
to Reissner-Nordstrom background, for which the analytic continuation of the metric into 
Euclidean space is given by 

ds'^ = f{r)dT^ + f-\r)dr^ + r^dO^ + sin^ ed<^'^ , (3) 

where / = {r — r-^){r — r^)/r'^, and using mass M and charge Q parameters of the black 
hole, we have 
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r± = M ± JM2 - Q2 , (4) 



For massless scalar fields the analytical approximation denoted by < (p >t reduces to 



K 



Therefore, in nearly extreme Reissner-Nordstrom spacetime such that /tr+ = {r^ — 
r_)/(2r_|_) ^ 1, the vacuum polarization of massless fields is strongly suppressed. (This 
is also justified by the result of Frolov estimated at the event horizon r = r+.) 

Such a excitation of vacuum polarization induced by finite black-hole temperature is an 
important aspect of quantum matter fields in black hole backgrounds, and it may remain 
valid for massive scalar fields too. Then, field mass m will just play a role of suppressing 
the amplitude of < 0^ > in comparison with massless fields. In this paper, however, we 
would like to emphasize another remarkable effect due to field mass, which we call mass- 
induced excitation as a remaining part of < 0^ > in the low-temperature limit T 0. Note 
that massive fields can have a characteristic correlation scale corresponding to the Compton 
wavelength 1/m. Our purpose is to show that nearly extreme (low-temperature) black holes 
can enhance the excitation of quantum fields with the Compton wavelength 1/m of order 
of the black hole radius (i.e., mr_|_ ~ 1). This mass- induced excitation may be expected as 
a result of wave modes in resonance with the potential barrier surrounding a black hole, 
for which the tail part of < 0^ > in the large-mass limit mr^ ^ 1 is generated with 
the amplitude decreasing in proportion to 1/iv? |^,|^ according to the DeWitt-Schwinger 
approximation developed by Christensen [Q. 

In this paper our investigation is focused on Reissner-Nordstrom background as the 
simplest example which allows us to consider the low-temperature limit Kr+ ^ 1 keeping an 
arbitrary value of mr^. (The black hole temperature and the field mass are measured in unit 
of the inverse of a fixed black hole radius r+. In Schwarzschild background with Kr+ = 1/2 
we cannot discuss the field-mass dependence of < 0^ > in such a low-temperature limit, 
and any resonance behavior of the polarization amplitude < 0^ > at mr_|_ ~ 1 will become 
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obscure by virtue of a contamination of the temperature- induced excitation in the mass range 
of mr+ ^1 [0].) Then, following the analysis given by Anderson and his collaborators 0,0, 
we compute the vacuum polarization of massive scalar fields, for which we have the analytical 
approximation of the form 

< 0^ >ap = < 0^ >T + < 0^ >m2 , (6) 

Here the additional contribution from field mass becomes 

27- (7) 

with Euler's constant 7. Unfortunately, this field-mass term contains a logarithmic diver- 
gence at the event horizon r = r_|_. Therefore, in Sec. II we develop the technique of 
analytical calculation to cancel such a divergent term, by paying the price that < cf? > 
is evaluated only near the event horizon. It is checked in Sec. Ill that the renormalized 
value of < 0^ > at the event horizon becomes identical, up to the leading terms of or- 
der of with the result derived by DeWitt-Schwinger expansion in the large-mass 
limit. In Sec. IV, using the small-mass approximation mr+ <^ 1, we show the tendency of 
temperature-induced excitation to be suppressed with incerase of field mass. We find in Sec. 
V the mass- induced enhancement of the polarization amplitude < 0^ >, by giving explic- 
itly the dependence on field mass in the low-temperature limit Kr+ <^ 1. The final section 
summarizes the results representing a remarkable difference of field-mass dependence of the 
polarization amplitude for scalar fields in equilibrium at various black-hole temperatures. 



II. CORRECTION TO THE WKB APPROXIMATION 

Let us start from a brief introduction of the method to compute the renormalized value 
of < 0^ > in Reissner-Nordstrom background (0), which has been developed by Anderson 
and his collaborators |0J^. Using Eq. (|I]) for a massive scalar field obeying the equation 

(□ - m^)(t){x) = , (8) 



the unrenormalized expression is given by 

< (fij-) > = lim{— ■ V c„cos(n/te)A„(r)} , (9) 

where Cq = 1/2 and c„ = 1 for n > 1. The separation of two points in Ge{x,x') is chosen 
to be only in time as e = r — r', and the radial part An{r) for each quantum number n is 
given by the sum of radial modes Pniij) and qni{r), 

OO -I 

Mr) = E{(2^ + ^)Pniir)qniir) - -y=} , (10) 
1=0 r\J J 

where / is the angular- momentum quantum number, and the subtraction term \jr^ is 

necessary for removing the divergence in the sum over /. The radial mode q^i satisfies the 

equation 

d^qr,i ^ 1 d{r^f)dqr,i ^ u^k"^ ^ /(/ + l)+mV% ^^ ^ ^ 
dr"^ r^/ dr dr /r^ " ' 

and it is chosen to be regular at r = cx) and divergent at r = r+. The same equation is 
satisfied by pnu which is chosen to be well-behaved at r = r+ and divergent at r = oo. 

The WKB approximation for the modes may be used to calculate the mode sums (p!OD, 
by assuming the forms 



and 



W = (2.W)i/2 ^^P(/(W/)^0, (12) 



(2rW)V2 

where the zeroth-order solution is chosen to be 



^ni = 777^4^ exp(- /(H^//)rfr) , (13) 



~ n\2 + {(/ + l)2 + ^V}4 • (14) 



To renormalize < 0^ > in the limit e — of point splitting, we subtract the counterterms 
< (p"^ >DS generated from the DeWitt-Schwinger expansion of < 0^ >, 



where o is equal to one half the square of the distance between the two points x and x', and 
0"" = V^a. Then, for the renormalized value defined by 

< 0^ >ren = < 0^ > _ < 02 ^ ^^^^ 

we can arrive at the analytical approximation (j^), if the second-order WKB approximation 
for W is used in the mode sums for n > 1 [0,0. 

Though Eq. (|]) can clearly show a spatial distribution of the vacuum polarization, the 
validity is rather restrictive. For example, in the asymptotically fiat region r — > oo it fails 
to give the expected dependence on field mass. It is instructive for later discussions to 
calculate precisely < 0^ >ren of thermal fields in equilibrium at a temperature T in fiat 
background (corresponding to / = 1), following the method of the Euclidean space Green's 
function Ge{x, x'). Denoting T by k/2ti, we obtain the exact solutions for p^i and Qni in flat 
background as follows, 

Pni = -i/j+i {rV m? + n^K^) , (17) 



and 



qm = ^Ki^{rVm^ + n^K^) , (18) 



and the mode sum over / in A„ results in 



An = -Vm2 + n2/t2 . (19) 

If we use the Plana sum formula for a function g{k) 
°° 1 dx 

T.9ij) = ^9ik) + g{x)dx + ij^ ^^^^ _ ^ \g{k + ix) - g{k - ix)] , (20) 



i=k 

the unrenormalized value is written by the integral form 

f"°o 2dn 

The first term in the right-hand side of Eq. (^) is completely canceled by the subtraction 
of the DeWitt-Schwinger counterterms ([T5|), in which we have a = — while the second 



< (P > = limj — -[— / dn co?,{nKe)\/m? + n^K^ + ( — V k^tt? — m?W . (21) 



term gives the renormalized value < 0^ >ren in flat background, which for massless flelds 
reduces to 

< 0^ >ren = TV12 , (22) 

and becomes equal to Eq. estimated in the asymptotically flat region. However, in the 
large-mass limit m ^ k, we obtain 

< 0^ >ren = m'l\T jll^f^-^l^ , (23) 

because the second integral over n in Eq. (|2l|) should run from the large lower limit mj k'^ \ 
to inflnity. This leads to a crucial difference from the approximated form (^, for which 
is expressed in inverse powers of uk such that 

^ -y + (^-m2)/2nK, (24) 

as a result of the mode sum over I using the zeroth-order solution (|T^ for W . It is clear 
that the sum of such an expansion form of An over n > 1 misses the exponential behavior 
g-27rm/K ^ 02 >^^^ jn asymptotically flat region. 

Now let us turn our attention to vacuum polarization at the event horizon / = 0, which 
is the main concern in this paper. Fortunately, we can claim that the above-mentioned 
deviation of Eq. (^ from the precise estimation becomes irrelevant, if we consider the limit 



f —> 0. This is because owing to the redshift factor f m W the expansion (|2J) remains 
valid even for a large mass m > k, by keeping the condition ra^f] j k <^ 1. Then, concerning 
vacuum polarization of massive flelds at the event horizon, we can use Eq. to show the 
dependence of < 0^ >ren on m. Of course, one may point out another crucial problem that 
Eq. contains a logarithmic divergence at r = r^. However, this singular behavior is due 
to the sum of An over the limited range of n > 1. Because the expansion form ( p^ also 
breaks down for n = 0, the contribution of Aq to < 0^ >ren is omitted in the calculation 
of Eq. (^. We would like to clarify an important role of the n = mode to give a regular 
value at the event horizon for the renormalized vacuum polarization < 0^ >ren (and also for 
the renormalized stress-energy tensor < Tat >ren)- 



To this end we propose the procedure to treat more precisely the mode sum over / in An 
at the event horizon, which is apphcable to the lower n modes. Note that near the event 
horizon the exact solution for should have the expansion form 
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Qni = In z asz' + E f^sz' , (25) 

with some coefficients as and jSg- The rescaled radial coordinate z is defined by z = (r — 
r^)/r^ ^ 1- This expansion form is not useful to calculate An at the event horizon, 
because the sums over / should be done without expanding in powers of z for requiring 
the convergence. Then, the important points to be mentioned here are the existence of the 
logarithmic term 2;"/^ In z and the power-law behavior z""/^ dominant for n > 1 in the limit 
z —>■ (except for the n = mode in which the logarithmic term becomes dominant). For 
the modes pni regular at the event horizon the dominant power-law behavior is given by 
z"/^, and the WKB forms ([13D and ([12D for Qni and Pni remain exact up to these dominant 



power-law terms. Hence, the value of An for n > 1 is exactly given by the WKB calculation 
in the limit z —>■ 0, and we will obtain a precise value of < 0^ >ren at the event horizon by 
taking account of the additional correction Aq to Eq. (|^). 

To resolve the problem of logarithmic divergence, however, it is important to note that 
the WKB form for fails to give the logarithmic behavior, which should play the role of 
canceling the logarithmic term contained in the DeWitt-Schwinger renormalization coun- 
terterms. (Because the leading logarithmic behavior in An would be z"' In z, the value of 
< 0^ >ren cau becomc regular at the event horizon only by considering a more precise treat- 
ment of the n = mode beyond the WKB level, while the same analysis for the n = 1 mode 
is also necessary to obtain a regular value of < >ren-) Hence, our key approach is to 
study the modified Bessel forms for the modes instead of the WKB forms as follows, 

Pni = {^f'Ux) , (26) 

and 

qm = {^f'Knix) , (27) 



where we have 

X = I {w/f)dr , (28) 



for which it is easy to check the vahdity of the Wronskian condition 

dqni dpni 1 

Pni-, Ini—r- = — ■ (29) 

ar ar r'^j 

The ordinary WKB forms are given if we assume Pni and g„/ to be proportional to /1/2 and 
respectively. Now, the function w introduced in place of W should satisfy the equation 



— {1 + — n^--} = — + ^ i 
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^ 1 d'^w 3 1 ^dwy ^ 1 d{r'^w) df 
2w dr^ Aw"^ dr 2r'^fw dr dr 



If w is rewritten into 



w = f/'y/r+ , (31) 
the solution of Eq. (^) allows the expansion form 

y = B{l + J2y,z') . (32) 

s=l 

From the well-known behavior of the modified Bessel function Kn{x) near x = 0, it is easy 
to see that qni has the expected logarithmic behavior near the event horizon. 

By substituting Eq. ( ^21) into Eq. (|30|) with the expansion in powers of z, we obtain the 
recurrence relation between the coefficients B and ys- For example, the lowest relation leads 
to 

^(n^ -l)(y,-2 + -^) = u{u + 1) + 2Kr+ - , (33) 

where z/(z/ + 1) = /(/ + 1) + m^r^. From the expansion up to the next power of z the relation 
between yi and 1/2 turns out to be 

2hir 



{n' - 4)1/2 = + l)z/i - l{l + 1) + f/(/tr+, n, y^) , (34) 
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where t/ is a slightly complicated quadratic function of yi which depends on n and nr^ only. 
An important point of the expansion form (0) is that we can require ys to remain finite 
in the limit / ^ oo, for which from Eqs. ( p3D and (Q) the asymptotic values of B and yi 
reduce to 

= l{l + 1) + m'rl + \ + n\2Kr+ -\) + 0(1-^) , (35) 

and 

yi = -i + o(r^) , (36) 

This dependence of ys on / allows us to calculate the mode sum over / in An by neglecting 
the terms with the higher powers of z in Eq. (^), and in the following Eq. ( p5D will be 
verified in terms of the cancellation of the logarithmic divergence in < 0^ >ren- 

We also remark that the amplitude of < 0^ >ren at the event horizon should not be 
interpreted as a quantity determined only by local geometry. The relations ( |33D and ( |3^ 
allow us to give a conjecture that the recurrence relation is truncated within a finite sequnce, 
and for the n-th mode the finite set consisted of B, yi, ■ ■ ■, yn-i is completely determined 
for any value of /. However, the coefficient yn remains unknown, unless the higher infinite 
sequnce of the recurrence relation is consistently solved for satisfying the boundary condition 
y — > (m^r^ + ra^K^r^)^/^ at 2; — > cxd as an eigenvalue problem. In particular, for n = we 
cannot give B for lower values of / without a further analysis of Eq. ([TT|) . This is the problem 
to be solved in the subsequent sections, and in this section we use Eq. ( pSj) for n = to 
derive the logarithmic term in Aq. 

By taking the limit z — 0, we can give the mode sum over / for n = written by the 
form 

°° 2/ + 1 I I 1 

= j:{^KoiB^2z/KrMB^/2z/Kr^) - (37) 

Then, we apply the Plana sum formula (^Dj) to Eq. (P7|), in which the modified Bessel 
functions is allowed to reduce to 
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Ko{Bpz/Kr+) ~ -7 - \n{B^z/2t,r+) , (38) 

and 



Io{BJ2z/Kr+) ^ 1 , (39) 



except for the integral defined by 



roo 2 + 1 I / 1 

/„ '''{^/fo(Bv^)^o(S\/^) - ■ («) 

To calculate the integral (^), let us recall that i? is a function of / satisfying 



2BdB/dl = 2/ + l + 0(r^) (41) 

in the large I limit and replace the integral of the modified Bessel functions over I by that 
over B to use the integral formula 

J 2BKo{Bv)Io{Bv)dB = B'^{Ko{Bv)Io{Bv) + Ki{Bv)Ii{Bv)} (42) 

for any variable v. Then, the same approximations with Eqs. (|38|) and ( p9D is applicable to 
the remaining integral given by 

dl , , „dB, 
and we arrive at the final result for Aq in the limit z ^ such that 



dl dR I I 

/ -^(2/ + 1 - 2B—)KoiB^2z/Kr^)IoiB^2z/Kr+) , (43) 



Ao = 4 + -{7 + ^ln(^)}, (44) 
K,r^ K 2 2Kr+ 



where 



1 r°o dR 

So = {Bl - )\nBo-^- dl{2l + 1 - 25—) \nB 
Z Z Jo al 

- ^^^m + l)\r.B{tl) + {2tl-l)\nB{-tl)} , (45) 
Jo e^'^' — 1 

if we denote B{1 = 0) by Bq. Hence, by adding kAq/Stt'^ to < 0^ >ap, the logarithmic 
divergence at the event horizon turns out to be canceled, and we obtain the renormalized 
value denoted by < 0^ >h as follows. 
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<<I>'>H= + ^{1 - In(mV^))} + . (46) 

It is interesting to note that the absence of the logarithmic divergence of < 0^ >ren at the 
event horizon is assured only by giving the asymptotic value (^) of B for the n = mode 
with very large /, which is determined through the local analysis near r = r+. Though in 
general we cannot obtain the renormalized value itself without deriving B for lower / modes, 
the large-mass limit can be an exceptional case for which the local analysis remains useful, 
and we calculate < 0^ >h up to the order of in the next section as a simple application 
of the procedure presented here. 

III. THE LARGE-MASS LIMIT 

To calculate the integral of B in 5*0 over / under the large- mass limit mr_|_ 3> 1, it is 
convenient to give the expansion form of B in inverse powers of + 1), by keeping the 
quantity fi = m^r^/z/(z/ + 1) to be of order of unity. (For the first integral present in 5*0 
we cannot assume Z(/ + 1) to be much smaller than mr+, while for the second integral the 
approximation yU ~ 1 — /(/ + l)(mr+)~^ may be allowed.) The expansion of B^ should be 
done up to the terms of order of l/z/(z/+ 1) for obtaining the terms of < cfP' >h- Then, 
the recurrence relation subsequent to Eqs. (|33D and (|3^) becomes necessary, for which the 
leading terms turn out to be 

2/2 = -| + ^(l-/x) + 0(m-2). (47) 

The key point of Eq. (^) is the absence of in the leading-order relation, from which Eqs. 
(B3) and (piD for n = can give 



= -1 + ^ + ^j^^^r, + 0(m-^) , (48) 

and 

= + + 2Kr+^) + ^^^^^V + 0{m-^) , (49) 
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where 

1 ,4 1 , 37 2 / X 

Now it is easy to calculate the integrals in Eq. up to the terms of order of (mr+)~^, 
and we can confirm the cancellation of all the terms much larger than (mr+)~^ in the 
expression (^Bf) for < (j? >h, giving the result 

<<P'>H= ^^r, I 2 4 (16^'^+ - 4/^^+ + 1) , (51) 

Note that the well-known term < cfp' >m-2 of the DeWitt-Schwinger approximation for 
< 0^ > can be written by 

< 02 > _ R^bcd _ R :[iah\ (52) 

for Reissner-Nordstrom background (with vanishing Ricci scalar), where Rabcd and Rab are 
the Riemann and Ricci tensors, respectively. If evaluated at the event horizon r = r+, this 
DeWitt-Schwinger term is found to be identical with Eq. (plj). Hence, for very massive fields 
with mr_(_ ^ 1 in equilibrium at black-hole temperature T = njl-R ^ we can claim the validity 
of the DeWitt-Schwinger approximation near the event horizon, as was previously shown in 
numerical calculations . Further, if mr+ is fixed, the tail part (0) in the range mr+ ^ 1 
becomes minimum at the black-hole temperature corresponding to ktj^ = 1/8, rather than 
at the low-temperature limit Kr+ ^ 1. The m-K coupling can give a slightly complicated 
change to the amplitude of vacuum polarization. In the next section we see a result of the 
m-K coupling as the suppression of temperature-induced excitation in a small-mass range. 



IV. THE SMALL-MASS LIMIT 

Now we consider scalar fields with very small mass mr+ <C 1, for which the temperature- 
induced excitaion given by Eq. (|^) will dominate. To reveal some correction due to the small 
field mass, let us begin with a brief analysis of purely massless fields. It is easy to see that 
Eq. ([TlD for the massless n = modes becomes equal to Legendre's differential equation, if 
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we use the variable x defined by x = 1 + (z/nr^). Then, from the behavior of Legendre 
functions at a; — 1 and x — > cxd, the modes qoi and poi should be chosen to be 

qoi = Qi{x) , poi = Pi{x) , (53) 

The mode sum in Eq. (0) for n = is known to be precisely zero for any x [§, and from 
Eq. (^ID the vacuum polarization at the event horizon reduces to 

which should be interpreted to be purely induced by the black-hole temperature. For purpose 
of extending the result to massive fields, it is useful to check explicitly through the procedure 
given in the previous sections that Sq in Eq. (|46|) vanishes. 

Recall that the function Qi{x) has logarithmic branch point at x = 1, and the dominant 
behavior near the point is 

Qi ^ lln(-l^)-^(l + 0-7, (55) 

where ip{s) is the logarithmic derivative of the gamma function (i.e., a polygamma function), 
and we have = —7 for Euler's constant 7. By comparing the logarithmic behavior of 
Qi with Eq. ( ^8]) for the modified Bessel function, we can determine the coefficient B as 
follows, 

B = exp{^(l + /)} . (56) 

To calculate the integrals over I in 5*0, we use integral representations for the polygamma 
function. For example, we obtain 

-{{2il + l)i){l + il) + {2il-l)i){l-il)] = 



e^'"' - 1 



/■°° , ,e-* 2t-2 + t-i l,cosh(t/2) , , 



4 sinh''(t/2) 
by virtue of the formula 
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Another useful formula is given by 

through which we arrive at the result 

^;|2e2^(i+oMi±i) _(2/ + i)}^^(i + /) = 
dl 

- + 7 6-2^--+/ dt(- + h7 + - • 60) 

Then, it becomes easy to calculate the integral over t for the sum of Eqs. (0) and (]60|), and 
we obtain Sq = 0. 

For the massive n = mode we rewrite Eq. (^) into the form 

{x^ - 1)^ + 2x^ - {/(/ + 1) + mV^ («:r+x + 1 - Kr+^jgoi = , (61) 

which can clarify the deviation from Legendre's differential equation. In this section a small- 
mass field having mr^ ^ 1 is assumed, and the solution perturbed by the field mass is given 
by 

qoi = Qi'ix) + qi{x) , (62) 

where / = 5 = 0{m?r']^). Because the terms proportional to m?r'^ in Eq. (|6TD is dependent 
on x, we use the recurrence formula valid for Qii (and also for P^) such that 

(/' + l)Qv+i - {21' + l)xQv + I'Qv-i = , (63) 

and the perturbed part qi is expanded in terms of Legendre functions as follows, 

Qi = Y.icfQv+k + c%Qv^k) ■ (64) 

k=l 

The coefficients Ck and C-k together with the eigenvalue 5 are determined by solving the 
recurrence relation 

cf{{l' + k){r + k + l)~ l{l + 1) - m^rlvfl) = m\l E(^f/J,cf|,. + 4724^) , (65) 
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where Cn^ = 1, and 



(0) n2 , 2 2 2i(2i + 1) - 1 

f, = 1 — Kr+ + K r,- , 

^ +^ +(2i- l)(2i + 3) ' 

vf^ = 2Kr+(l - Kr+)^T^ , v^i'^'^ = 2/tr+(l - /tr+)^-^ , (66) 

(2) 2 2 + + 2) (_2) 2 2 



(2^ + 3)(2z + 5) ' ' +(2i-3)(2i- 1) ■ 

Then, the first-order perturbation is found to be 

m'^nrl Kr+ {I + + 2)Qi+2 l{l - l)Qi-2 ^^ 



and 



5 = — ^{ 1 - Kr+ 2 + kV — ^ — -J- -} , (68) 

2/ + l^^ +(2/ - 1)(2/ + 3)^ ' ^ ^ 



for which the coefficient B is estimated to be 



ciZ + +^ /(/ + !) + (2/-l)(2/ + 3)^^ ' ^^^^ 
Using these equations, one may calculate the polarization amplitude < 0^ >h at the event 
horizon. However, for I = the value of B becomes divergent as a result of the existence 
of the undefined function Q_k in Eq. (0). This will mean a dominant contribution of the 
/ = mode in the small-mass limit. 

To estimate more precisely B = Bq for I = 0, the subscript / in the Legendre functions 
should be replaced by taking account of the approximate relation Qs-k — Pk-i/^ for 
5^1. Then, the term m^r'^Qs^i which appears in go should be interpreted to be of order 
of unity, contradictory to the perturbation scheme. This problem is resolved if we add 
another independent solution for Eq. (|6T|) written by 

oo 

Po = S^^Ps + Y.{df^Ps+k + d^-lPs-k) (70) 



fc=i 



to go as follows. 



Qo = J2i(^k^Qs+k + C-lQ5-k) +Po , (71) 

k=l 
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where we require that 5 ^c?| + dg'^ = £ <S 1 for dg'^ of order of unity. Of course, the 
coefficients df^ should satisfy the same recurrence relation with c^j^\ and we obtain for 
k > 1 

4ti = Oiimr^r) , = Oiimr^r) , (72) 

in addition to the ratio d^^l/d^j^l^ = 0{m'^r'^). Then, the asymptotic behavior of the / = 
mode Qoo at x ^ 1 is approximately given by 

1 ^ r(/c + (1/2)) 1 (0) ,(0)W2 Xfc (r,o\ 

which should be consistent with the boundary condition 

X 

goo — — exp(— m^ria;) (74) 

X 

at a distant region far from the event horizon. 

To check the consistency, let us derive the approximate recurrence relation which is valid 
up to the leading order of m^r^ and reduces to 

r^°^ d^°'^ 2 A- - 1 

~ 4t2 " ^ (2^ + mk - mk - 3) ' ^''^ 

and 

S-'c% + dfl, - '^^'H2k + 2){4k + l){4k-l)- ^''^ 
Noting the relations between the lowest coefficients such that 

5-icL°l = 2m'^Krl{l- Kr+) (77) 



and 



-iJO) _L ^(0) _ ^2 ^2^4 



+ d^ = , (78) 



we arrive at the result 
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which can satisfy the boundary condition if e = —rriKr'^. 

Unfortunately, we cannot determine e to the order of m^r^, unless the reccurence relation 
is studied to the higher order. Hence, we only keep the leading correction of order of mr_|_ 
in the I = mode, 

qoo - Qo- mnrl , (80) 

which means that Bq = e^'^(l + rriKr^). For the I > 1 modes Qqi we must also consider the 
perturbation with the terms written by the Legendre functions Pk{x). However, it is sure 
that no perturbation of order of mr+ does not appear for / > 1, and we obtain 

S'o — —ln(l + mKri) ~ — m^ri , (81) 



if we omit the higher-order corrections. Now the vacuum polarization given by Eq. ( PB| ) for 
small-mass fields becomes approximately 

<<P'>H^ ^^(l-3mr+), (82) 

which cleary shows that the temperature-induced excitation is suppressed by field mass. 
As m becomes larger, the amplitude may monotoneously decrease in the whole mass range 
extending to mr^ ^ 1 where the DeWitt-Schwinger approximation < 0^ >h^ {mr^)~'^ is 
valid. This simple dependence on m is supported through numerical calculations for several 
values of mr_|_ in Schwarzschild background (/tr+ = 1/2) ||^. In the next section, however, 
we point out a different dependence on field mass, which is a resonant behavior of < 0^ >h 
remarkable in the low-temperature case Kr+ ^ 1. 



V. MASS-INDUCED EXCITATION 

Let us turn attention to quantum fields at the event horizon of nearly extreme black 
holes to show an interesting feature of the mass-induced excitation of vacuum polarization. 
Then, we do not limit the range of the parameter mr+, but we solve Eq. (^Tj) under the 
assumption Kr+ ^ 1 by the help of the technique of asymptotic matching. 

18 



At large values of x Eq. (^) reduces to the form 

-fT- + — 3 ( 2 + +m K r^)qoi = , (83) 

in which we cannot neglect the terms depending on Kr+ to require the exponential decrease 
of goi- For the approximate differential equation we obtain the solution 

go; = W_^^^^^^^{2mKr\x)lx , (84) 

where Waf, denotes the Whittaker function with the asymtotic behavior 

WaAu) ~ M'^exp(-M/2) (85) 

as M — »• 00. This asymptotic soluition can remain valid in the range 

1 < X < l/fi;r+ , (86) 

where we obtain the approximate behavior 

goz ^ ^^"^'^ ~ ^\ {2mKr\xY^^x-^ + ^^^^ ^ ilmnrXxy x'^ . (87) 



Note that if a; <^ Eq- ( PH ) becomes approximately equal to Legendre's differential 

equation, giving the solution 

goi = CP,[x) + Dg,(x) . (88) 

The coefficients C and D should be determined by the matching with the approximate 
solution (pTl), and it is easy to see that the ratio C jD is of order of (mKr^)^*^"*"^. Hence, we 
can neglect the term in goi, and the asymtotic behavior at x — > 1 turns out to be 

go^ ^ -D{^ln(^)+7 + ^(z/ + l)} , (89) 

from which we obtain 

B = e'^(^+^) , (90) 

for calculating Sq (and < 0^ >h) through Eq. 
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A useful expression of 5*0 to understand the field-mass dependence is derived if we use 
the integral formula 

+ 1) = 2 + " + 4' + io (e^.< + l)(t^ + .'^ + . + (l/4)) • P'l) 

In fact, for F{1) = (—i){(2il + 1) lni?(i/) + {2il — 1) In B(—il)} which is one of the integrands 
in 5*0, we obtain 

,2 ^^2 , ,21 , . , ^ ^ r 8^^^ /2 + (l/2)-t2-( 



= /ln{(/^ - C)^ + I'} + arctan(^) - ^2.^ + 1 _ ,2 _ ^)2 + p ' (92) 

where ( = m^r^ + (1/4), and the value of arctan(u) runs from to vr in the range < u < 00. 
Further, the integral given by 

J dl{2l + l-2B^)\nB (93) 

is rewritten into the form 

^-Mu + 1) + \}{Hu{u + 1) + ^) - 1} - e^^('^+i){V^(z. + 1) - ^} 



+2( ^^^Ht' + ^i^ + l) + -^), (94) 



which is equal to zero as / — > 00. We therefore arrive at the result 



where 



l.> Is, > C r°^r tG(t) Hit) , , 



2 , 1 ^ dl P + (1/2) - - ( 



and 



iy(t) = tln{(t2-C)2 + t2|^^^^^^^(^_^ _ (97) 

Under the low-temperature approximation Kr+ ^ 1 we neglect the term fi;/247r^r+ in Eq. 
(HO), and the polarization amplitude at the event horizon is finally given by 



,22 i2 w?r\ 



In(^) - \{l + InC) + Hi-^nK + -^!rK}dt ■ (98) 

m^ri 8 Jo e^^^* + 1 e^^^* - 1 
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Now it is easy to check the value of < 0^ in the large- mass limit mr+ ^ 1, and we 
obtain 

for which we can reconfirm that it is equal to the DeWitt-Schwinger approximation (with 
/tr+ — > 0). We can also consider the small-mass limit mr+ <^ 1 under the condition m/zt ^ 1, 
and the approximate expression of < 0^ >h becomes 

Sn^rl < (f)^ >H ^ -mV^{^ + 7 + ln(mr+)} , (100) 

which can remain positive by virtue of the existence of the logarithmic term — m^r^ ln(mr+). 

We evaluate numerically the integrals in the expression of < 0^ >h, and the field- 
mass dependence is shown in Fig. |l[ Note that the maximum excitation of < 0^ >h 
occurs at mr+ ^ 0.38, and the peak amplitude denoted by < 0^ >max is estimated to be 
Svr^r^ < 0^ >^ax— 0.0424. We can clearly see a resonance behavior of the polarization 
amplitude for massive fields with the Compton wavelength 1/m of order of r+ and also the 
tail part given by Eq. (pQ]) in the mass range of mr+ ^ 1. 



VI. SUMMARY 

We have studied vacuum polarization of quantized scalar fields in Reissner-Nordstrom 
background by means of the Euclidean space Green's function. In particular, the renormal- 
ized expression < 0^ >h at the event horizon r = r+ has been derived by revealing the 
contribution of the n = mode, which can cancel the logarithmic divergence. 

We have found the dependence of < 0^ >h on field mass m: (1) The tail part observed in 
the large-mass limit mr+ ^ 1 becomes equal to the DeWitt-Schwinger approximation. (2) 
For small-mass fields a suppression of temperature-induced excitation due to the coupling 
between m and k, occurs according to < 0^ >h=< 0^ >r (1 — 3mr+), where the massless 
part with the amplitude proportional to the black- hole temperature T = k/27t is given by 
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Svr^r^ < 0^ >T= Kr+/3. We can expect that mass-induced excitation becomes important 
for massive fields with mr+ ~ 1. Unfortunately, it is difficult to investigate in detail various 
aspects of the m-n coupling in the case that both mr+ and kt^ are of order of unity. (3) Our 
main result therefore has been to show a resonance behavior of mass-induced excitation of 
vacuum polarization around nearly extreme Reissner- Nordstrom black holes with nr^ ^ 1: 
If the Compton wavelength 1/m of a massive field is of order of the black-hole radius r+, the 
amplitude of vacuum polarization has a peak at the resonance mass given by mr+ ~ 0.38. 

There should be a critical temperature Tc = Kc/2tt of black holes in the range < 
fi;r+ < 1/2, below which a resonance peak of < cfP' >h is observed in the field-mass de- 
pendence. (If K > Kc, the polarization amplitude monotoneously decreases with increase of 
m.) Though the value of Kc remains uncertain within the analysis presented here, it is sure 
that dominant fields as quantum perturbations near the Schwarzschild horizon should be 
massless, while nearly extreme holes will have a quantum atmosphere dominated by fields 
with a resonance mass. The peak amplitude given by Svr^r^ < 0^ >^cix— 0.0424 at the 
nearly extreme Reissner-Nordstrom horizon is not so smaller than the massless part given 
by Svr^r^ < 0^ >t= 1/6 at the Schwarzschild horizon with the same area Anr'^. (If com- 
pared under the same black-hole mass M, the former becomes slightly larger than the latter 
evaluated by Stt^M^ < 0^ >t= 1/24.) Considering a black hole evolving toward the zero- 
temperature state with a fixed radius r+, we conclude that the mass m of dominant fields 
generating vacuum polarization shifts from mr_|_ ^ 1 to mr_|_ ~ 0.38 as the contribution of 
mass- induced excitation becomes important, without changing the polarization amplitude 
so much. Quantum back-reaction due to massive fields will become very important for 
nearly extreme (low-temperature) black holes. 
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FIG. 1. The field-mass dependence of vacuum polarization < >h at the nearly extreme 

Reissner-Nordstrom horizon r = r+. The amplitude has a resonance peak at mr+ 2± 0.38 and a 

tail part decreasing in proportion to for very massive fields. 
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